Fugacity and Reheating of Primordial Neutrinos 
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We clarify in a quantitative way the impact that distinct chemical Tc and kinetic Tk freeze-out 
temperatures have on the reduction of the neutrino fugacity below equilibrium, i.e. Ti/ < 1, 
and the increase of the neutrino temperature Tv via partial reheating. We establish the connection 
between T^, and Tk via the modified reheating relation T^, (T„) /Tj , where T-y is the temperature of the 
background radiation. Our results demonstrate that one must introduce the chemical nonequilibrium 
parameter, i.e., the fugacity, T^, as an additional standard cosmological model parameter in the 
evaluation of CMB fluctuations as its value allows measurement of Tfc. 

PACS numbers: 51.10.-|-y,95.30.Cq,14.60.Pq,26.35.-|-c 



Introduction: The free-streaming relic neutrino 
distribution is an important input into structure for- 
mation in the universe and the calculation of CMB 
fluctuations Recent results from the nine-year 
WMAP observations. Table 7 of Ref.[3 and another 
independent study of BBN 0] favor an effective num- 
ber of neutrinos at BBN of = 3.55tJ^;4^ and 
= 3.71^045, respectively, while the newly re- 
leased Planck data finds ^ 3.30 ±0.27 H. How- 
ever, this fit produces a 2.5 s.d. tension with direct 
astrophysical measurements of the Hubble constant. 
Including priors from SN surveys removes this ten- 
sion and result in = 3.62 ± 0.25. 

Much work has been done to compute the modi- 
fication of the neutrino distribution due to two-body 
interactions between neutrinos and e^e~ in general 
relativistic kinetic theory. Standard electroweak in- 
teractions give rise to Ny = 3.046 0-01 while a study 
of modified neutrino-electron interactions leads to 
Ni, = 3.12 Q. Current observational data does not 
have tight enough error bounds to either confirm 
or show incompleteness of the standard picture of 
neutrino freeze-out, but the deviation of the cen- 
tral value determined by Planck and others from 



the computed value suggests that one consider ad- 
ditional factors capable of modifying the neutrino 
distribution and thereby increasing iV^. 

One such mechanism is the interaction of neutri- 
nos with collective degrees of freedom in the e^e~ 
plasma. The calculations referenced above assume 
that the freeze-out process is dominated by two- 
body scattering, but scattering against collective 
degrees of freedom is often a significant factor in 
plasma dynamics. If such an effect is relevant here, 
the result would be a lower neutrino kinetic freeze- 
out temperature and, as we will show, a larger 
value of . We do not attempt to model the multi- 
body scattering processes that could lead to a reduc- 
tion in Tfe. Rather, we perform a parametric study 
of the effect that the reduction of Tk has on the 
neutrino distribution and photon-neutrino reheat- 
ing ratio within a simple model that emphasizes the 
physics of interest. We discover a new effect that ac- 
cumulates during the temporary reduction of the de- 
celeration parameter q caused by the electron mass, 
which has a significant effect on the neutrino distri- 
bution when Tfc w m^. 

Our study relies on the following two proper- 
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ties of the dynamical evolution of the primordial 
6+, e~, u, 9, 7 plasma: 

1) The chemical and kinetic, often called ther- 
mal, freeze-out conditions are distinctively differ- 
ent. Once the universe temperature drops below the 
chemical freeze-out temperature Tc, there are no re- 
actions that, in a noteworthy fashion, can change 
the neutrino abundance. Wc will find that the pre- 
cise value of Tc is immaterial to the present study as 
long as > 2me , which is the consensus of all eval- 
uations [3, Il0| - in our own detailed study, we find 
Tc ~ 2.7me [lil. As the universe cools further, at 
some point neutrinos begin to free-stream [l2[ . This 
is the kinetic freeze-out and Tk is the key parameter 
in determining the neutrino momentum distribution. 
Note that particle creation, which determines is 
a hard process while exchange of momentum, which 
determines Tk, is soft. Hence Tc does not have the 
same complications involving collective phenomena 
as Tk does. This allows us to treat Tc as being de- 
termined by two-body interactions while Tk will be 
treated as an unknown parameter. 

2) There exists a brief time window when the tem- 
perature of the universe T < 2me, during which 
the deceleration parameter q drops mildly from that 
for the pure radiative universe because the electron- 
positron mass becomes significant and affects the dy- 
namics of the expansion. This lasts until the density 
of e'^e~ pairs becomes negligible and the fully radia- 
tive expansion resumes at T < O.lSmg. 

We find that these effects combine to produce 
a photon- neutrino temperature ratio T^JT^ that is 
closer to 1, a known 'neutrino reheating' effect [B-0|- 
However, it has thus far escaped attention that re- 
heating of neutrinos can be accompanied by a signif- 
icant underpopulation of neutrino phase space rela- 
tive to an equilibrium distribution, characterized by 
a little-known cosmological model parameter which 
will be called neutrino fugacity, Ti,. Its significance 
for neutrino cosmology has been previously recog- 
nized but is not widely appreciated. 

We relate T^jTy to T,y and discuss how a signifi- 
cant deviation of < 1 alters the neutrino momen- 
tum distribution at the time of recombination. Our 
findings will demonstrate that should be included 
among the standard cosmological model parameters 
governing the experimental CMB observables, espe- 
cially with regard to neutrino mass constraints. 

Neutrino Energy Spectrum and Moments: 

Consider a particle with degeneracy gp, mass m, and 
chemical freeze-out temperature Tc- Prior to chemi- 
cal freeze-out the distribution function has the usual 
Bose-Einstein or Fermi-Dirac chemical equilibrium 
form. After chemical but prior to kinetic freeze- 



out at Tk, collisions between particles maintain a 
maximum entropy distribution, but particle number 
changing processes are no longer able to occur. The 
number of particle-antiparticle pairs is therefore con- 
served. To achieve particle conservation in an envi- 
ronment in which temperature changes, another pa- 
rameter is needed: T. The momentum distribution 
and thus value of T is found by maximizing entropy 
at fixed particle number, antiparticle number, and 
energy. 

The thermal (kinetic) equilibrium momentum 
distribution functions arc applicable down to kinetic 
freeze-out Tk 

f ^9v_ 1 ^^ 

where T = 1 for T > and T 7^ 1 for T e [T^, Tk] 
is the same for both particles and antiparticles. One 
can also write 

^ = TlnT + ^t. (2) 

T controls the number of particle-antiparticle pairs 
available to fill the phase space. In the absence of 
(measurable) true chemical potentials /it originated 
from conserved quantities such as lepton number, T 
is identical to the complete particle fugacity. There- 
fore, in this situation an effective chemical potential 
H is derivable from T through Eq. ([2]) alone. Chem- 
ical non-equilibrium is characterized by T 7^ 1, i.e. 
for fit — but fJ, 0- For fermions the limit T ^ 1 
approaches a degenerate Fermi distribution, while 
for T ^ the Boltzmann limit may be used. 

The moments of Eq. ^ define the energy den- 
sity p, pressure P, number density n, and entropy 
density s of a distribution in kinetic equilibrium 
Eq. ([1]). Standard thermodynamic identities, such 
as the Gibbs-Durham relations, remain valid down 
to the kinetic freeze-out stage in the evolution of 
the universe. In principle, once the neutrinos start 
free-streaming these relations do not apply anymore. 
However, in a radiation dominant universe a signifi- 
cant violation occurs only when the mass of the neu- 
trino becomes a significant scale. 

Modeling Universe Expansion: The neutrino os- 
cillations are fast compared to the dynamics of the 
universe for T ~ T^. ~ mg. Thus each neutrino 
spends one-third of its time in each of the three 
flavor states and its kinetic reactivity will be the 
arithmetic average of the three relevant rates. Addi- 
tionally, neutrinos and antineutrinos in a practically 
matter-antimatter symmetric universe are subject 
to the same dynamics. Finally, at the temperature 
scale of interest the neutrino masses are negligible. 
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Together, these factors imply a common value of 
for all three flavors of neutrinos and their antipar- 
ticles. Hence we can collect them all into a single 
distribution function, with a degeneracy gi, = 6. 

The three relevant dynamical quantities are thus 
the temperature T of the 7, e+e~, v, V system, the 
common neutrino fugacity T^, and the scale factor 
a of a homogeneous, isotropic, and spatially flat uni- 
verse: 

= - a{tf{dx'^ + dy"^ + dz'^). (3) 

The dynamics are determined by: 

a) The Einstein equation. For a flat FRW universe 

we integrate 



■ \ 2 

a 



3M2 ' 



Mp = l/VS^, (4) 



where p is the total energy density of photons, e~^e~ , 
and neutrinos. 

b) Divergence freedom of the total stress-energy ten- 
sor (which is implicit in the Einstein equation) . Us- 
ing the Gibbs-Durham relations, we find that diver- 
gence freedom of the stress-energy tensor. 



v^r^o = p + 3(p + p)- = o. 



is equivalent to 



(5) 



(6) 



c) A dynamical equation governing the neutrino- 
number-changing process. 



which reads 



(7) 



(8) 



Here — >■ 1. For further details of the derivation 
of Eq. dil) see Eq.(7.290) in Ref.Q, of which this is 
a slight modification, as well as Refs.fisl. [T6l|. 

Neutrino oscillations allow us to collect all three 
neutrino flavors into a single distribution for which 
the dominant reaction strength is that in Eq. ([7]), 
which is the reaction most responsible for chem- 
ical equilibrium as it is the only rate that con- 
tains charged current (exchange of VF^), with small 
contributions by the neutral current (Z°) mediated 
e^e~ annihilation, which also feeds into v^^ and Vt- 



Since each of the three neutrino flavors can be fed 
in e^e~ annihilation the total relaxation time is 



f 



1 

Tt 



I 



(9) 



where the electron rate dominates within the preci- 
sion required in our approach. 

The rate constant Te(T) is obtained by con- 
ventional kinetic methods described in Ref . |14| : see 
Ref.[llj for more details. Our electroweak ampli- 
tudes and invariant rates agree with those presented 
in Rcfs.jl, We obtain, to better than 5% preci- 
sion in Boltzmann approximation, the result 



1 -f- (2me/T)2e0-29(2m,/T3Tg)- 



\ (10) 



where Gp is the Fermi coupling constant. The neu- 
trino degeneracy in the last normalizing term in 
Eq. ([TU]) is f/i, = 2 X 3/ . Note that the flavor fac- 
tor 3 / averages the rate constant shown in Eq. (|9]) 
and thus the effective rate of chemical equilibration 
is 1/3 of that for i/e alone in absence of neutrino 
oscillations. 

This kinetic theory model describes a smooth 
chemical freeze-out process, as is inherent in Eq. ([5]). 
This is in contrast to the instantaneous transition 
approximation, in which at some temperature Tc it 
is assumed that all neutrino number changing pro- 
cesses cease completely. Under this 'sudden' approx- 
imation, the two terms in Eq. ^ are separately zero 
and and Eq. (|6]) and Eq. ^ are replaced by the inde- 
pendent conservation of comoving entropy and neu- 
trino number. 

In the more precise 'smooth' model of Eq. ([5]), 
comoving entropy is not exactly conserved. In gen- 
eral, particle number changing processes occurring 
outside chemical equilibrium, i.e. T ^ 1, are entropy 
generating processes. We verified that in our case 
in ultrarelativistic plasma, the conversion of €^e~ 
pairs into neutrino pairs is practically entropy con- 
serving. This explains why under the entropy con- 
serving approximation, T^(T) remains accurate to 
within a few percent, a point we will demonstrate 
explicitly. Thus, in principle there is no need to 
model the chemical freeze-out in detail for the pur- 
pose of the discussion carried out below and we are 
justified in assuming exact conservation of entropy 
when convenient. This would of course be anyone's 
first guess but considering the issues of reheating we 
address, we have convinced ourselves by modeling 
the smooth chemical freeze-out that this is an accu- 
rate model assumption. 
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FIG. 1: Neutrino fugacity as a function of neutrino 
kinetic freeze-out temperature Tk/nie- The solid line 
corresponds to the exact numerical solution for and 
the dot-dashed line closely following it is the entropy 
conserving freeze-out approximation; see the text for de- 
tails. Dashed line (scale on right) is the deceleration 
parameter. 

Effect of Electron Mass: If all relevant particles 
have negligible mass then both entropy and number 
density scale with T^, hence the solutions to Eq. ^ 
and Eq. ([8]) are simply 

T, = 1, Toe i. (11) 
a 

However, the proximity of Tc and Tk to the electron 
mass modifies the situation. As the temperature 
passes through the electron mass, the deceleration 
parameter in the expansion of the Universe, 



displays a small deviation from the pure radiation 
dominant value of q = 1, decreasing slightly towards 
the matter dominant value of q = 1/2, indicating 
a temporary period of matter relevance. It is this 
drop in q that causes T^, to decrease and the neutrino 
phase space to become underpopulated. Though the 
magnitude of the drop in q (right scale in figure [T]) 
is not large, the effect on is cumulative. 

We obtain a quantitative measure of this under- 
population effect by numerically solving Eq. dH]), 
and ^ to obtain T,y{Tk) for each choice of the neu- 
trino kinetic freeze-out temperature Tk, which we 
treat as a free parameter. We begin at a tempera- 
ture high enough, e.g., T = lOme, for the universe 
to be fully radiation dominant and neutrinos to be 
in chemical equilibrium with — 1, and integrate 



across Tc until T = Tk, at which point the neu- 
trinos become free-streaming. The result is shown 
as a solid line in figure [TJ We also show the re- 
sult of the entropy conserving chemical freeze-out 
approximation as a dot-dashed line, which closely 
parallels the solid line. In the entropy conserving 
approximation, the calculation is insensitive to the 
precise value used for the neutrino chemical freeze- 
out temperature T^, as long as it is several times 
the electron mass, i.e., larger than the temperature 
where begins to deviate from 1. This is so since 
neutrino chemical freeze-out occurs in the nearly ex- 
act radiation dominant era of the universe and prior 
to the temperature where the influence of the elec- 
tron mass on the universe dynamics takes effect. In 
this domain the neutrino distribution evolves as in 
Eq. irrespective of whether or not it is in chem- 
ical equilibrium. The difference between these two 
calculations is less than 1%, so for the purposes of 
this work, the approximate model is valid. 

Delaying the thermal freeze-out condition until 
it is close to or below the electron mass results in 
neutrinos acquiring energy liberated by e~^e~' anni- 
hilation via scattering and thus are reheated, along 
with photons, down to Tk . Once thermal freeze- 
out occurs, only photons are reheated and from that 
point on T^ > Ty. At Tk neutrinos begin to free 
stream and the neutrino entropy is conserved in- 
dependently of the entropy in photons and e^e~, 
as shown in Ref.[l2|. We make the approximation 
that the total comoving entropy is exactly conserved, 
hence whatever entropy remains in e"*"e~ at Tk goes 
solely into reheating the photon temperature. 

The calculation of the resulting photon-neutrino 
temperature ratio T^/T^ is the same as that used to 

1/3 

derive the traditional reheating ratio of (11/4) ' , 
only that here one starts with less comoving entropy 
in e^e~ since by the time T = Tk, some of the en- 
tropy has already been transferred to neutrinos. We 
carry out this procedure numerically and find a sim- 
ple power-law relation to better than 1% between 
the reheating ratio T^/T^ and Tj^: 

/ rp \ 2.57 

= 0.420 (^^ j . (13) 

Note that for ^ 1, we obtain the standard re- 
heating ratio Try/Tjj — (11/4)^/^, as expected. Re- 
ducing Tk leads to a reduction of T^jTy compared 
to standard reheating, and introduces a significant 
deviation of Ti, from unity, an effect which has not 
been recognized before. 

Momentum distribution at recombination: In 
Ref.[l2| we consider and solve the Einstein- Vlasov 
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equation for the free-streaming form of the neutrino 
momentum distribution after kinetic freeze-out. Us- 
ing the fact that nii, <C 71-, the result is 



Mt,p) - 



1 



87r3 T^^eP/^- + 1' 



(14) 



where the neutrino temperature is obtained by 
solving Eq. p3p . This is a nonequilibrium distribu- 
tion for two reasons: 

a) The appearance of an effective chemical potential 
Eq.©: =T^lnT^. 

b) The absence of the neutrino mass in the exponen- 
tial means that when the temperature approaches 
and drops below m,y, the distribution is no longer of 
the kinetic equilibrium form. 

Depending on the ratio of neutrino mass to re- 
combination temperature, one or both of these ef- 
fects may be significant. From Eq. (|T4|) we see that 
the number density is calculated as if the neutrino 
were massless, but to compute the energy density or 
momentum one weights the energy or momentum of 
a massive particle with an effectively massless dis- 
tribution function at a temperature T^. A fugacity 
T,y 7^ 1 together with the modified reheating ratio 
/T^{T u) imply that the effect of the neutrino mass 
being close to the recombination temperature is al- 
tered from the standard result when is reduced. 



Interpretation of Planck Results To ex- 

plain a value N^, = 3.62 ± 0.25 at ion-electron 
recombination, we find Tk 



O.lQtoM MeV and 



= 0.71_g gg. The value of effective chemical po- 
tential /iy at recombination, Eq. ([2]), follows, noting 
that is preserved in the free-streaming neutrino 
momentum distribution and hence it continues to 
be present at the time of recombination. We find 
//^ ~ -0.086 eV. 

The above estimates explicitly assume that the 
effect we present here is the only one available to 
explain the Planck results (with priors) for 7V^. Fur- 
thermore, we reinterpret the fitted Ni, rather than 
performing a fit in which ^ \ and the reheat- 
ing constraint Eq. p3l) is included. The neutrino 
momentum distribution that would enter such a fit 
would be a background gas with a smaller neutrino 
number density and greater momentum than cur- 
rently assumed. Incorporating together with the 
modified reheating ratio, Eq. ([T3|) . into fits of CMB 
data has the potential to constrain the magnitude of 
the neutrino mass with greater consistency. 

To interpret the literature value = 3.046, we 
require r^(°) = 0.806 MeV, which leads to an effective 
chemical potential = T^lnT^ ~ —2.2 x 10~^T^. 
The value Tk we find to interpret a value Ni, ~ 3 is 



very sensitive to the method of achieving the kinetic 
decoupling; it corresponds to T^°^ = 0.7864 MeV 
obtained by computing directly the two-body col- 
lision freeze-out condition [ll|. The theoretical un- 
derstanding of how a reduction from the lowest order 
standard model two-body scattering value T^°' ~ 0.8 
MeV to Tk ~ 0.2 MeV is possible, is a question re- 
quiring further study. 

Summary and Conclusions: We have derived the 
dependence of the background cosmic neutrino dis- 
tribution and reheating ratio on the neutrino ther- 
mal freeze-out temperature Tk- The drop in the de- 
celeration parameter of the universe when the elec- 
tron mass becomes a relevant scale, T ~ nie, com- 
bined with distinct neutrino freeze-out temperatures 
that satisfy the condition Tc > me > Tk, where Tc 
is the neutrino chemical freeze-out temperature, re- 
sult in a deviation of the neutrino distribution from 
equilibrium. This effect is described by a primordial 
neutrino fugacity different from unity. 

The standard cosmological model determines 
neutrino freeze-out by two-body microscopic scat- 
tering 0, 01 . However, currently a value of 3.6 > 
Ni, > 3 is favored by BBN related observations 0,0 
and the recent Planck fit Q with priors. The pro- 
posed modification of Tk from a two-body scatter- 
ing value is motivated by the mechanism of neutrino 
scattering against collective degrees of freedom in 
the e"'"e~ plasma, but in this work Tk was left as a 
free-parameter . 

We determined the magnitude of the fugacity 

as a function of Tk, sec figure [TJ Wc established 
a modified reheating relation between the current 
photon and neutrino temperatures required by T^, ^ 
1, Eq. (1121), and thus demonstrated a possible large 
modification of the primordial neutrino momentum 
distribution. should thus serve as a cosmological 
model parameter, hitherto presumed to be close to 
unity, in the evaluation of CMB fiuctuations. 

A value well above Ni, ~ 3 would support a 
lower kinetic freeze-out temperature Tk and hence a 
larger effective chemical potential in absolute value, 
up to ^i, ~ — 0.87Ti,, maybe required. The neu- 
trino distribution evolves after freeze-out as a free- 
streaming gas and this preserves the value of T^, 
implying that neutrinos have a temperature depen- 
dent effective chemical potential fXi, = T^ In Ti, after 
freeze-out. One must note that at the recombination 
temperature Tr = 0.25eV, this effective chemical po- 
tential is on the order of the neutrino mass limit and 
is therefore non-negligible. 
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